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EXECUTIVE  SUMMARY 


OBJECTIVE 

This  report  summarizes  the  basic  and  Brownian  Reber  search  theories  which  have 
been  used  to  model  the  performance  of  the  Advanced  Unmanned  Search  System.  It  is 
intended  to  serve  as  a  reference  to  the  system  modeling  which  influenced  the  design  of 
the  system.  It  should  also  serve  as  a  foundation  for  the  design  of  optimal  searches. 

APPROACH 

Reber’s  original  theory  has  been  modified  to  incorporate  navigation  and  control  fea¬ 
tures  typical  of  autonomous  vehicle  systems;  the  natural  modification  involves  modeling 
vehicle  trajectories  as  fractional  Brownian  paths.  The  general  structure  of  the  Brownian 
Reber  theory  is  identical  to  the  structure  of  the  basic  Reber  theory,  so  that  the  usual 
figures  of  merit,  such  as  mean  time  to  detection,  are  calculable  in  terms  of  system 
measurables,  such  as  navigation  error  and  lateral  range  function. 

RESULTS 

The  basic  and  Brownian  theories  are  described  in  some  detail,  and  their  similarities 
and  differences  are  explicitly  analyzed.  The  form  of  the  expression  for  mean  time  to 
detection  is  provided.  The  important  parameters  in  the  Brownian  theory  are  determined 
numerically  for  a  significant  range  of  system  variables,  and  these  results  are  presented 
in  a  series  of  figures. 
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INTRODUCTION 


The  Advanced  Unmanned  Search  System  (AUSS)  was  developed  by  the  Naval 
Command,  Control  and  Ocean  Surveillance  Center  (NCCOSC)  to  improve  the  Navy’s 
ability  to  find  and  identify  items  lost  or  placed  on  the  sea  floor  at  depths  as  grert  as 
20,000  feet.  Items  such  as  the  Palomares  H-Bomb,  the  U.S.S  Scorpion ,  the  U.SI-'. 
Thresher,  Korean  Airlines  Eight  007,  Air  India  Flight  182,  and  the  cargo  do.  r  of 
United  Airlines  Flight  811  are  examples  of  equipment  lost  by  the  U.S.  and  othu 
countries.  Searching  for  these  items  proved  difficult  and  highlighted  a  critical 
technology  area:  deep  ocean  search. 

The  first  purpose  of  a  search  theory  model  is  to  quantify  the  performance  of  a 
search  system  in  terms  of  the  parameters  which  most  directly  affect  that  performance, 
and  as  with  any  complicated  system,  this  quantification  necessitates  a  number  of 
important  simplifications.  The  object  of  this  report  is  to  explicitly  state  these  simplifi¬ 
cations,  derive  a  theory  based  on  these  simplifications,  and  provide  a  parameterized 
summary  of  the  performance  of  the  AUSS  which  may  be  used  in  its  engineering  and 
operation. 

In  the  continuing  engineering  of  AUSS,  the  acceptability  of  such  simplifications  is  a 
subject  of  constant  discussion,  particularly  when  the  optimal  operation  of  the  system  is 
being  considered.  This  report  is  intended  to  document  refinements  to  Reber’s  theory 
(Reber,  1956),  which  led  to  the  initial  design  and  development  of  AUSS.  These  refine¬ 
ments  are  suggested  by  a  better  understanding  of  AUSS  as  it  actually  exists,  as 
opposed  to  how  it  might  have  been  generally  designed,  and  while  it  is  accepted  that 
the  theory  which  is  the  subject  of  this  report  may  itself  be  refined  in  time,  it  is  com¬ 
plete  within  the  limits  of  the  simplified  search  scenario  outlined  in  the  next  section. 

The  next  section  of  this  report  outlines  the  important  parametric  simplifications,  the 
resulting  theory,  and  its  application  to  the  quantification  of  system  performance.  This 
report  then  discusses  the  AUSS  navigation/control  model  and  an  explicit  derivation  of 
the  location  probability  expressions.  The  RMS  path  deviation  and  mean  search  path 
length  are  also  derived.  Reber’s  theory  is  discussed  in  some  detail,  so  tint  it  will  be 
clear  where  it  differs  from  the  modified  theory  which  is  the  subject  of  1  iis  report. 
These  differences  are  discussed  with  the  object  of  substantiating  the  need  for  a  modi¬ 
fied  theory.  The  next  section  concentrates  on  determining  the  factors  which  influence 
the  time  it  takes  to  complete  the  nth  coverage  of  the  nominal  search  area,  which  in 
turn  enters  into  the  mean-time-to-detection  figure-of*merit.  Explicit  expressions  for 
various  search  scenarios  of  interest  to  AUSS  are  derived  and  explored.  The  final 
section  summarizes  the  influence  of  the  various  system  parameters  through  a  series  of 
figures. 


1 


SEARCH  THEORY  PARAMETERS 


SEARCH  SCENARIO 

For  AUSS,  the  simplified  search  scenario  is  as  follows:  The  search  vehicle,  or 
search-sensor  platform,  is  assumed  to  be  searching  a  rectangular  area  by  following  a 
sequence  of  parallel  paths  and  employing  a  search  sensor,  such  as  a  sonar  or  a 
camera,  which  periodically  provides  information  about  or  “searches"  rectangular  sub- 
areas  along  the  path.  In  this  description,  the  overall  performance  of  the  search  system 
is  the  result  of  its  performance  in  two  essentially  separable  activities,  that  is,  how  well 
the  system  follows  the  sequence  of  parallel  paths  (the  “nominal”  search  pattern)  and 
how  well  the  search  sensor  searches  the  subareas  to  which  it  is  periodically  applied. 
How  well  the  system  follows  its  nominal  paths  is  a  result  of  the  accuracy  and  precision 
of  the  navigation  data  it  uses,  as  well  as  the  control  algorithms  which  actuate  the 
vehicle’s  thrusters,  rudders,  elevators,  and  other  parts;  how  well  a  search  sensor 
searches  an  area  is  described  by  the  conditional  probability  of  an  object  of  interest 
being  detected  by  the  sensor,  given  that  it  is  known  to  be  present  in  the  sensor’s  field 
of  view. 

This  simplified  search  scenario  is  schematically  depicted  in  figure  1.  The  search 
vehicle  attempts  to  follow  the  nominal  paths,  of  length  L,  in  the  directions  indicated  by 
the  arrows  superimposed  on  the  paths;  and  at  regular  intervals  separated  by  a  distance 
8  along  the  nominal  paths,  a  search  sensor  is  employed  to  provide  information  about  a 
subarea  assumed  to  be  centered  on  the  position  of  the  vehicle.  These  subareas  have  a 
width  W  (across  track,  or  perpendicular  to  the  nominal  path)  and  a  length  8  (along 
track,  or  parallel  to  the  nominal  path),  the  dimension  W  being  related  to  the  range  of 
the  search  sensor.  The  nominal  paths  are  separated  by  a  distance  d,  so  that  the  nomi¬ 
nal  search  area  has  width  and  length  H  -  Md,  L  ■  N8,  respectively,  that  is,  there  are  M 
such  nominal  paths  and  N  subareas  along  each  nominal  path. 


2d  <M-1)d  x 

Figure  1,  Nominal  search  scenario. 
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NAVIGATION  ERROR 

The  model  assumed  to  describe  the  vehicle’s  success  in  following  the  nominal 
search  pattern  is  outlined  in  more  detail  below,  but  in  short  it  states  that  the 
probability  that  the  perpendicular  distance  of  the  vehicle  from  the  path  will  change 
from  one  seprch  interval  to  the  next  by  an  amount  e  is  given  by 


Note  that  the  second  moment  of  this  distribution,  o,  is  a  function  of  the  distance 
between  search  swaths,  8,  along  the  nominal  path,  although  in  the  following  expres¬ 
sion,  this  dependence  will  usually  be  suppressed.  In  any  event,  this  second  moment  is 
the  so-called  navigation  error  of  the  system.  Most  importantly,  this  model  results  in  the 
square  root  of  the  mean  perpendicular  deviation  from  the  nominal  path  (the  RMS 
deviation)  being  an  increasing  function  of  distance  along  the  path,  with  the  specific 
form 

**•«>  -  (£),/J» 

where  y  is  the  distance  along  a  nominal  path  (0  <  y  <  L).  Note  that  the  parameters 

{v,6,y} 

are  measureables,  so  that  this  expression  gives  the  search  theory  parameter  o  as  a 
function  of  the  demonstrated  capabilities  of  an  actual  system. 

LATERAL  RANGE  FUNCTION 

In  general,  the  probability  that  a  search  sensor  will  detect  an  object  of  interest 
when  it  is  known  to  be  in  the  field  of  view  of  the  sensor  is  a  function  of  the  object,  its 
location  in  the  field  of  view,  the  nature  of  the  environment  in  which  the  sensor  is 
being  employed,  and  the  performance  of  the  operator  or  algorithm  making  use  of  the 
sensor  information.  The  first  simplification  is  to  restrict  the  field  of  view  of  the  sensor 
by  indicating  its  effective  range,  which  in  the  present  situation  is  a  perpendicular  dis¬ 
tance  (the  lateral  range)  from  the  sensor  beyond  which  it  has  no  chance  of  detecting 
an  object  of  interest.  For  AUSS,  the  search  swath  is  assumed  to  be  symmetrical  about 
the  position  of  the  sensor,  so  that  the  area  searched  by  the  sensor  in  one  application 
has  width  W  and  length  8,  as  above,  where  the  sensor  has  lateral  range  W/2.  The  sec¬ 
ond  simplification  lumps  the  remaining  factors  into  a  single  number,  a  conditional 
probability  (3,  which  describes  the  probability  that  a  given  object  in  a  given  environ¬ 
ment  will  be  detected  by  a  given  operator  or  algorithm  provided  that  it  is  within  the 
field  of  view  of  the  sensor.  This  is  the  so-called  “cookie-cutter"  lateral  range  function, 
and  though  quite  crude  it  is  a  useful  rendition  of  the  actual  sensor  performance.  Later 
reports  will  consider  the  effect  on  overall  search  performance  of  allowing  p  to  be  a 
true  function  of  lateral  range,  as  in  p  ■  P(X),  where  X  is  the  lateral  range  of  the  object 
of  interest  from  the  sensor. 
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LOCATION  PROBABILITY 


The  probability  that  an  object  of  interest  will  be  detected  after  one  attempted 
search  of  the  nominal  search  area  by  following  the  sequence  of  parallel  paths  is  the 
location  probability,  <p .  There  is  a  canonical  form  for  <p  which  can  be  developed  as 
follows.  Suppose  that  the  search  vehicle  follows  the  nominal  search  pattern  with  no 
error  and  detects  objects  in  the  field  of  view  of  its  sensor  with  probability  3,  as  above, 
and  suppose  further  that  the  object  of  interest  is  certain  to  be  located  at  least  some* 
where  within  the  nominal  search  area.  Since  the  search  proceeds  with  no  navigation  or 
control  error,  it  effectively  partitions  the  nominal  search  area  into  NM  subareas  of  size 
W8.  Now,  the  probability  that  the  object  of  interest  is  located  within  any  one  of  these 
subareas  is  the  ratio  of  the  magnitude  of  the  subarea  to  the  magnitude  of  the  whole 
search  area,  namely,  W8///L,  so  the  probability  of  not  detecting  the  object  on  one 
nominal  coverage  of  the  search  area  is 

l  “  hl) 

Since  the  location  probability  is  the  probability  that  the  object  will  be  detected  after 
one  nominal  coverage 


or 

* "  i  -  «p(-^)  - 1  -  «p(-^r) 

since  8  -  L/N  and  d  *  HIM.  It  is  shown  in  the  next  section  of  this  report  that  within 
certain  reasonable  limits  on  the  relative  values  of  the  parameters  discussed  thus  far, 
the  location  probability  for  a  system  which  does  not  search  with  perfect  navigation  and 
control  is  well-approximated  by 


where 


K  - 


The  forms  of  the  location  probability  and  the  other  parameters  are  dictated  by  the 
desire  to  be  consistent  with  Reber’s  theory.  The  parameter  \0  must  be  determined 
numerically,  but  in  practice  it  is  a  convenient  measure  of  the  effect  of  imperfect  navi¬ 
gation  on  Overall  search  performance  as  quantified  by  the  location  probability  0 . 


4 


MEAN  TIME  TO  DETECTION 

A  real  search  is  made  up  of  a  number  of  activities  whose  execution  contributes  to 
the  overall  tin.e  a  search  is  said  to  require;  these  activities  include  those  things  which 
may  be  cor.-sidered  to  directly  contribute  to  the  probability  of  finding  an  object  of  inter¬ 
est  as  well  as  those  things  which  are  more  logistic  in  nature.  In  the  former  class  is  the 
search  activity  itself,  that  is,  the  use  of  the  search  sensor  to  provide  information,  while 
in  the  latter  category  are  such  things  as  transit  time  to  the  search  area  and  system 
refurbishment  time,  in  which  it  is  assumed  that  the  system  has  a  limited  deployment 
time  which  is  less  than  the  time  required  to  conduct  a  single  search  of  the  nominal 
search  area.  Similarly,  usually  the  case  that  search  first  proceeds  with  a  sensor  which 
views  large  areas  rather  quickly,  such  as  a  sonar,  but  which  cannot  distinguish  between 
the  object  of  interest  and  a  set  of  other  objects  with  similar  signatures  (false  targets). 

In  this  case,  a  contact  evaluation  sensor  is  employed  to  distinguish  between  true  and 
false  contacts.  The  use  of  this  second  sensor  is  known  as  contact  evaluation,  and  the 
time  spent  on  contact  evaluation  depends  upon  the  density  of  false  targets  in  the 
nominal  search  area  as  well  as  the  time  it  takes  to  evaluate  a  given  contact,  whether 
the  contacts  are  evaluated  as  they  occur  or  in  groups  at  some  later  time. 

It  is  important  to  note  that  it  is  assumed  that  more  than  one  nominal  coverage  of 
the  search  area  will  be  attempted  before  an  object  of  interest  is  located.  Denote  by 
T(n)  the  time  spent  on  all  activities  during  the  nth  attempted  coverage  of  the  nominal 
search  area.  Observe  that  in  terms  of  the  location  probability  <f>,  the  probability  that 
the  object  of  interest  will  not  be  located  until  the  nth  attempted  coverage  is  given  by 

m  -  (i  -  0-v 

so  the  mean  time  spent  searching  for  an  object  of  interest  is 

t-  i mm  -  i n»xi  - 

n  ■  1  n  ■  1 

This  mean  time  to  detection  is  a  useful  figure  of  merit  through  which  the  effects 
may  be  investigated  of  the  various  system  parameters  on  search  system  perfomiance. 
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BROWNIAN  REBER  SEARCH  THEORY 


NAVIGATION  BY  DIRECTED  RANDOM  WALK 

Consider  a  nominal  search  scenario  as  outlined  above  and  represented  schemati¬ 
cally  in  figure  1.  The  AUSS  vehicle  attempts  to  follow  the  “square-wave”  nominal 
search  pattern  by  using  a  Doppler  velocity  sensor  whose  integrated  output  provides 
estimates  of  net  changes  in  position  in  the  search  area;  this  position  information  is 
used  by  Type  1  control  loop  algorithms  which  attempt  to  keep  the  vehicle  on  the  nomi¬ 
nal  path  at  every  point.  The  net  result  of  integrating  a  sensor’s  output  (the  Doppler 
velocity  estimates)  and  using  a  line-following  control  algorithm  based  on  integrated 
sensor  output  is  modeled  here  by  assuming  that  the  probability  of  the  perpendicular 
distance  of  the  vehicle  from  the  nominal  path,  evaluated  at  regular  intervals  along  the 
nominal  path,  changing  by  an  amount  e  from  one  interval  to  the  next,  is 

p(e)<,e  ■  sWexp('w)e 

where  8  is  the  separation  of  intervals  along  the  nominal  path. 

In  terms  of  figure  1,  in  which  a  position  along  the  nominal  path  is  described  by  the 
y-coordinate  and  a  perpendicular  position  relative  to  the  nominal  path  x  -  x0  is 
described  by  x  -  x0 


p{e)d€  m  p{xn\xn.x)dxn 


C*n 


Here  the  index  n  runs  from  n  *•  1  to  n  ■  N,  where  N  -  LI 8,  L  being  the  length  of 
the  nominal  paths,  which  are  successively  described  by  x0  -  qd ,  the  index  q  running 
from  qmQtoqmM-l,  where  M  -  Hid,  H  being  the  width  of  the  search  area  and  d 
the  separation  of  the  nominal  paths. 

The  probability  of  the  vehicle  following  a  path  described  by  a  specific  sequence  of 
perpendicular  positions  {X\,  Xi,  ....  x^}  at  the  intervals  {8,  28,  ....  JV8)  along  the 
nominal  path  is  then 


p{xx\x0)dxx  •  p{x2\xi)dx2 . p(xN\xN-i)dxN 


2if  J 


dx\dx2...dxN 


(1) 


Equation  1  is  the  fundamental  directed  random  walk  assumption  which  determines 
the  navigation/control  performance  of  the  system. 
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LOCATION  PROBABILITY 


Now  consider,  as  above,  a  search  sensor  with  a  cookie-cutter  lateral  range  function 
.vith  swath  width  W  and  detection  probability  (5.  If  the  object  of  interest  is  located  at 
(jc, ,  y,)  within  the  nominal  search  area,  then  y*  ^  y;  <  y*+i  for  some  \  <  k  <  N.  The 
probability  of  detecting  the  object  while  searching  the  nominal  path  x  =  x0  is  then 


PM 


(*n  " 


2cr 


since  the  object  is  detected  with  probability  (3  when  x,  -  W/2  s  x  £  y,  +  W/2  and 
y  ■  yk,  and  zero  otherwise.  Straightforward  integration  with  changes  of  order  and 
variable  gives 

w 

■  iki  I  w  'xpf-^h  (2) 

F'nally,  observe  that  the  probability  of  failing  to  detect  the  object  while  searching 
the  nominal  path  x  »  x0  is  1  -  Pk  (x0),  so  the  probability  of  failing  to  find  the  object 
on  one  nominal  coverage  of  the  search  area  is 


Equations  2  and  3  define  the  location  probability,  <f>(xt,  y,);  a  more  useful  quantity 
is  the  average  of  <p(xt,  y,)  over  all  possible  locations  of  the  object  within  the  nominal 
search  area.  Provided  that  the  magnitude  of  the  area  of  the  paths  along  the  edges  of 
the  search  area  is  small  relative  to  the  magnitude  of  the  total  area,  averaging  on  x, 
can  be  performed  over  the  smaller  range 


x'  -  d/2  s  xt  <  x'  +  d/2 

for  some  x'  well-removed  from  the  edges  of  the  nominal  search  area.  This  gives  the 
average  location  probability 


or,  more  explicitly 


<P 


1 

dL 


2 


j  j  <p(x„  y,)dx4y, 


(4) 


0  s 


1 

dL 


I 


(1  ~Pk(Xo  = 


7 


Now  define 


L  2  Af-1 


X  =  - 


pw 


r  \  t  r  1 

Xr\dL  f  J  n 

^  r\  rt  0  —  0  + 


(5) 


o  x,  £  =  o 

X  ’  2 


then  <f>  has  the  familiar  form 


PW , 


<t>  **  l-exp(-~-A)  ««  l-exp(-mA) 
a 

with  m  -  $W/d,  the  nominal  path  density. 

It  is  possible  to  considerably  simplify  Eq.  5  for  K;  first  write 
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Ihe  argument  of  the  exponential  in  the  latter  integrand  varies  slowly  over  the  range 
of  integration  for  a  considerable  number  of  the  likely  values  of  the  parameters  of 
interest  here;  it  is  shown  later  in  this  report  over  exactly  what  range  of  these 
parameters  this  argument  may  be  replaced  by  its  average  over  the  nominal  search 
area.  Replacing  the  argument  by  its  average  gives 
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leading  to  an  approximation  of  N. 
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and  an  approximation  of  <p 


<po  ■  1  -  exp(-mAa) 


It  remains  only  to  simplify  the  expression  for  X0.  Notice  first  that  the  inner  integral 
over  x,  is  constant  for  y*  £  yt  £  y*+i,  so  that  the  integral  over  y,  is  actually  a  sum: 
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each  term  in  the  sum  being  multiplied  by 

y*4i 

J  dy,  =  y*+i-y*  =  <5 
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Using  the  explicit  form  for  Pk  (jt0)  and  interchanging  the  sum  on  q  with  the  integra¬ 
tion  over  x  gives 
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Making  the  change  of  variable  z  -  x,  -  qd  for  each  term  in  the  sum  over  q  gives 
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Since  x  **  (H  -  d)/2  is  the  ^-coordinate  most  distant  from  the  edges  of  the  search 
area,  using  x1  ■  (//  -  d)/2  in  the  last  expression  gives  the  form  discussed  above: 
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RMS  PATH  DEVIATION  AND  MEAN  PATH  LENGTH 

The  square  root  of  the  mean  of  the  square  of  the  perpendicular  deviation  of  the 
actual  vehicle  path  from  the  nominal  path  is  a  function  of  the  distance  along  the 
nominal  path.  With  these  successive  distances  given  by  yk  »  £8,  as  above,  where  the 
index  k  runs  from  k  »  1  to  k  «  N,  the  RMS  path  deviation  is  found  from 

* 


(  1  v  7  7  r  -■ 

1  - 1  — J*- 


dxk 


9 


Set 


/ 


(72b)  i 


Explicit  evaluation  trivially  gives 
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Changing  the  order  of  integration 
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or,  more  loosely 
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Note  that  it  is  implicit  here  that  each  nominal  path  is  begun  correctly. 

In  the  parametric  studies  of  this  report,  the  RMS  deviation  will  be  summarily  repre¬ 
sented  by  the  constant 
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It  is  similarly  straightforward  to  find  the  mean  path  length,  with 
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giving  the  mean  length  of  a  single  nominal  path,  so  that  the  total  mean  path  length  for 
one  nominal  coverage  is  the  sum  of  M  such  quantities.  It  is  easy  to  see  that 
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where 
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It  is  possible  to  find  an  asymptotic  form  for  /,  which  is  reasonably  accurate  for 
8  «  L  and  o  <  8,  using  the  method  of  steepest  descents  (Jeffreys  and  Jeffreys,  1956): 


and 


The  approximations  of  10  can  be  used  to  further  illustrate  the  differences  between 
Reber’s  basic  theory  and  a  Brownian  Reber  theory.  In  the  Brownian  theory,  the  vehicle 
paths  are  well-known  fractals  (Mandelbrot,  1982),  and  the  fractal  dimension  of  these 
paths  serve  to  reveal  these  fundamental  differences.  The  expression  i0  is  the  length  of 
a  segment  of  the  vehicle’s  path,  so  the  finite,  D-dimensional  measure  of  an  entire  path 
of  N  m  L/8  such  segments  is 

Ld  m  lim  ■  lim 

i  -*  a  0  6-*  o  0 

where  the  D-dimensional  measure  of  each  segment  is  nD{ 8)  ■  /?( 8). 

The  two  approximations  of  J0  given  above  distinguish  two  fundamentally  different 
types  of  path.  For  o  <  8 
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Now 


so  if  D  >  1 
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By  presumption  Ld  is  finite,  soZ)>l-»o**<5  2 ,  which  together  with  o  <  <  8 
requires  D  £  1,  a  contradiction.  Hence,  for  o  <  <  8  as  8  -*  0,  the  path  has  dimension 
D  m  l  and  the  limiting  path  length  is  L. 

On  the  other  hand,  for  a  >  8 


and  a  similar  limiting  analysis  shows  that 
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lim  4o°l 

6-*  o  O 


Again  the  presumed  finitude  of  LD  requires  a  ~  6D  or  o  -  tjdD  for  some  r\  £  1, 
which  is  the  description  of  a  fractional  Brownian  path  (Mandelbrot  and  van  Ness, 
1968)  of  dimension  1  £  D  $  2.  The  parameters  T|  and  D  are  determined  from 


log  a  -  log  log  8 

using  linear  regression  on  the  measurables  o  and  8;  note  that  the  fractal  path  length  is 
•  then 

L,  . 

The  difference  between  the  two  types  of  path,  namely,  those  for  which  o  <  8  and 
o  >  8,  is  essentially  the  difference  between  those  paths  whose  pointwise  variance 
vanishes  faster  than  the  decreasing  scale  at  which  it  is  evaluated  and  those  whose 
pointwise  variance  vanishes  more  slowly.  The  latter  paths  are  the  model  chosen  here, 
since  choosing  the  former  implies  that  the  vehicle’s  path  evaluated  at  discrete  times  is 
smoother  than  in  practice  it  appears  to  be  (Walton,  1992a,  1992b). 

Two  fundamental  differences  between  the  basic  Reber  theory  and  a  Brownian  Reber 
theory  are  explicit  now,  the  first  already  identified  above:  the  RMS  path  deviation  is  an 
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increasing  function  of  distance  along  the  nominal  path,  which  is  characteristic  of  sys¬ 
tems  which  emplby  Type  1  control  loops.  In  the  basic  theory,  the  navigation  error  and 
the  RMS  path  deviation  are  identical  and  constant  along  the  path,  as  indicated  in  the 
summary  discussion  of  the  basic  theory  given  below.  The  basic  Reber  path  is  conse¬ 
quently  two-dimensional,  while  the  above  argument  shows  that  paths  in  a  Brownian 
theory  have  dimension  1  £  D  £  2,  where  D  is  determined  experimentally.  This  second 
fundamental  difference  between  the  basic  and  Brownian  theories,  as  well  as  the  need 
to  model  navigation  using  Type  1  control  loops  shown  by  the  first  difference,  led  to 
the  Brownian  theory  which  is  the  subject  of  this  report. 

Figure  2  shows  a  nominal  Reber  path  and  three  Brownian  Reber  paths  for  three 
fractal  dimensions  1  <  D  <  2;  all  four  paths  are  of  equal  probability.  Figure  3  shows 
the  path  multiplier 


as  8  -4  0  with  r\  •  1. 


NOMINAL  REBER 

Figure  2.  Equiprobable  nominal  Reber  and  Reber-Brownian  paths  of  fractal 
dimension  D. 


13 


i - 1 - 1 - 1 - 1 -  ■  i - r  u 

-3  -2.6  -2  -1.6  -1  -0.8  0 

Iog10(6/LENQTH) 

Figure  3.  Path  multiplier  as  a  function  of  o0  . 

REBER’S  THEORY  IN  OUTLINE 

To  a  great  extent,  Reber’s  theory  of  search  for  narrow  path  locators  underlies  all  of 
the  above,  and  since  this  report  is  intended  to  document  a  significant  modification  to 
this  basic  theory,  it  is  appropriate  here  to  present  it  in  outline.  The  main  points  of 
divergence  between  the  modified  and  basic  theories  will  be  noted  as  they  are 
developed. 

The  fundamental  difference,  from  which  all  others  follow,  lies  in  the  form  which 
the  navigation  error  is  assumed  to  take.  For  the  basic  theory,  the  probability  that  the 
search  vehicle  is  within  «  of  the  nominal  path  is  the  same  as  the  probability  in  the 
modified  theory  of  a  change  of  magnitude  c  in  the  x-eoordinates  of  succeeding  posi¬ 
tions.  As  the  vehicle  attempts  to  follow  the  nominal  path  x  ■  x0,  for  any  y-coordinate 
along  the  path,  the  x-coordinate  has  probability 

** "  (9) 

where  the  navigation  error  is  now  quantified  by  the  parameter  5. 

This  is  to  be  compared  with  Eq.  1.  Note  that  in  the  simplified  search  scenario  dis¬ 
cussed  above,  navigation  error  for  the  modified  theory  is  a  function  of  the  along-track 
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evaluation  interval  6.  In  the  basic  theory,  the  navigation  error  is  independent  of  any 
such  parameter,  and  as  a  consequence,  the  RMS  path  deviation  is  identical  to  this 
navigation  error.  Figure  2  shows  four  typical  paths,  of  identical  relative  probability,  for 
these  differing  basic  and  Brownian  navigation  assumptions.  It  is  the  demonstrated  simi¬ 
larity  of  the  path  for  the  modified  theory  to  actual  AUSS  paths  and  the  assumed  con¬ 
sequences  of  the  use  of  a  Doppler  velocity  sensor  whose  integrated  output  provides 
estimates  of  relative  position  for  line-following  Type  1  control  loop  algorithms,  which 
led  to  the  development  of  the  Brownian  Reber  Theory. 

For  a  cookie-cutter  lateral  range  function  with  a  conditional  probability  (3,  with  the 
identical  simplified  search  scenario  discussed  above,  the  probability  of  detecting  the 
object  while  searching  the  nominal  path  x  -  x0  is 

w 

P>M '  /  „  exp('i?)fa  (10) 

which  is  to  be  compared  with  Eq.  2.  This  leads  to  a  location  probability  in  the  equa¬ 
tion 


[1-P(x0-O)][l-P(x0md)]...[l-P(x0m  (M-l)d)]  ■  l-0(x„y,)  (11) 

which  may  be  compared  with  Eq.  3.  This  is  then  averaged  over  ^-coordinates  of  object 
locations  in  the  range  x'-d/2  £  xt  £  x'  +  d/2  to  give 

i-  ri  ii -p(*o-«<oi 

q  ■  0 

Notice  that  the  modified  theory  also  averages  the  location  probability  over  the  y- 
coordinate  of  object  locations;  in  the  basic  theory,  the  y-coordinate  plays  no  role. 

Now  defining 


the  location  probability  takes  the  canonical  form 

<p  a  1-exp^-^yj 

Equation  12  may  be  compared  with  Eq.  5.  The  simplification  of  the  expression  for 
y  proceeds  in  essentially  the  same  fashion  as  that  of  k  leading  to  the  approximation 
Yo' 
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This  last  expression  is  to  be  compared  with  Eq.  6. 

As  indicated  above,  die  RMS  path  deviation  for  the  basic  theory  is  identical  to  the 
navigation  error,  <5 .  Since  the  navigation  error  is  not  a  function  of  the  along-track 
evaluation  interval,  8,  a  mean  path  length  is  not  calculated  since  it  leads  to  unbounded 
estimates  of  mean  path  length  in  the  basic  theory. 

In  any  event,  once  a  location  probability  is  available,  the  mean  time  to  detection  is 
calculated  in  precisely  the  same  fashion.  The  utility  of  Reber’s  theory  is  due  in  part  to 
the  parametric  summarization  of  the  dependence  of  the  location  probability  on  naviga¬ 
tion  error,  path  separation,  swath  width,  and  lateral  range  function,  presented  in  a 
series  of  figures  which  accompany  the  documentation  of  the  theory  (Reber,  1956).  A 
similar  summarization  is  included  in  a  later  section  of  this  report. 
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FIGURE  OF  MERIT 
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TOTAL  TIME  FOR  THE  nth  COVERAGE 

As  discussed  earlier,  a  real  search  consists  of  a  number  of  activities  which  contrib¬ 
ute  to  the  total  time  a  search  is  said  to  take,  including  the  time  spent  transiting  to  the 
search  area,  time  spent  actually  employing  search  sensors  to  provide  information,  time 
spent  evaluating  contacts,  and  time  spent  on  support,  logistic,  and  general 
housekeeping  activities,  during  which  the  search  is  said  to  be  in  abeyance.  For  the 
moment,  the  time  spent  initializing  the  search  will  be  taken  for  granted  and  will  not 
appear  in  any  expressions  for  mean  time  to  detection.  Since  it  is  assumed  that  more 
than  one  coverage  of  the  nominal  search  area  will  probably  be  needed  before  an  object 
is  detected,  denote  by  T(n)  the  total  time  spent  on  the  nth  coverage;  then  summarizing 
these  remarks 

T(A0  -  Ts(n)  +  TC(V)  +  TA(n) 

where  T,  ( n )  is  the  time  spent  actually  employing  search  sensors  to  provide  informa¬ 
tion,  or  the  search  time,  Tc(n)  is  the  time  spent  investigating  contacts,  or  the  contact 
evaluation  time,  and  TA(n)  is  the  time  spent  on  all  other  activities,  or  the  abeyance 
time.  With  <P  the  location  probability,  the  mean  time  to  detection,  also  discussed 
earlier  is 

00  00  00 

T  -  £  T(n)  P(n)  .  2  r(/iW(l-0"-'  £  [7i(n)  +  Tc(«)  +  W)(l  -0)"” 

n  ■  1  « *  1  » ■  1 


SEARCH  TIME 


Assuming  that  the  nominal  search  pattern  consists  as  above  of  M  parallel  paths  of 
length  1  and  separated  by  a  distance  d  and  that  the  vehicle  while  actually  searching 
moves  with  speed  v,  the  time  spent  employing  search  sensors  to  provide  information  is 
proportional  to  the  number  of  coverages,  as  in 


Ts(n)  *  n 


Ml 


where  E  is  the  mean  path  length  discussed  earlier. 


CONTACT  EVALUATION  TIME 


It  is  assumed  that  the  search  proceeds  first  by  using  a  sensor  which  has  a  large 
field  of  view  (a  broad  area  search  sensor,  such  as  a  sonar)  but  which  cannot  distin¬ 
guish  between  the  object  of  interest  and  a  set.  of  objects  which  have  similar  signatures, 
i.e.,  false  targets;  these  contacts  are  evaluated  by  using  a  separate  sensor  (such  as  a 
high-resolution  CCD  camera).  The  total  time  spent  on  contact  evaluation  is  then  a 
function  of  the  number  of  false  targets  present  in  the  nominal  search  area,  quantified 
by  the  false  target  density,  /,  and  the  time  spent  investigating  each  contact,  quantified 
by  the  contact  evaluation  time,  t.  It  is  assumed  here  that  each  contact  is  immediately 
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evaluated  and  that  this  leads  to  its  elimination  as  a  potential  contact,  whether  true  or 
false.  This  assumption  is  known  as  immediate  contact  evaluation  with  elimination. 

Since  a  false  target  is  indistinguishable  from  the  object  of  interest  by  the  broad 
area  search  sensor,  it  has  the  same  location  probability,  <t>,  as  the  object  of  interest, 
and  since  the  nominal  search  area  has  magnitude  HLt  the  number  of  false  contacts 
generated  on  the  first  attempted  coverage  is  taken  to  be  (fHL) 4>.  The  number  remaining 
to  be  detected  on  the  second  coverage  is  then 

fHL  -  fHLtf)  -  fHL(\  -<f>) 

which  leads  to  fHL(\  -  <p)<t>  false  contacts  on  the  second  attempted  coverage.  In 
general,  then,  die  time  spent  evaluating  contacts  after  n  attempted  coverages  of  the 
nominal  search  area  is  taken  to  be 

tm  - 

r  ■  l 


ABEYANCE  TIME 

Since  a  typical  nominal  search  area  can  be  as  large  as  100  sq  nmi,  it  must  be 
assumed  that  the  search  vehicle  will  be  deployed  and  recovered  many  times  during  a 
nominal  coverage.  In  particular,  for  AUSS  a  certain  , amount  of  time  per  deployment  is 
spent  on  descending  to  search  depth,  ascending  from  search  depth,  and  maneuvering  at 
search  depth,  maneuvering  consisting  at  the  very  least  of  turning  from  one  nominal 
path  to  the  next.  Likewise,  each  deployment  is  likely  to  involve  a  refurbishment  of  the 
vehicle  (battery  changeouts,  for  example),  as  well  as  other  logistics.  Suppose,  for  each 
deployment,  that  Tx  is  the  time  spent  for  descent,  T2  the  time  spent  for  ascent,  TM  the 
time  spent  maneuvering  at  depth,  7*  the  time  spent  for  vehicle  refurbishment,  and  TL 
the  time  spent  on  other  logistics,  and  assume  that  these  quantities  remain  constant  for 
each  deployment.  If  No  is  the  number  of  deployments  in  any  one  nominal  coverage, 
assumed  to  remain  constant  for  each  coverage,  the  abeyance  time  is 

Un)  -  nND(T,  +  T2  +  TM  +  7>  +  TL) 

Note  that  all  these  quantities  are  specified  by  the  system,  with  the  exception  of  the 
maneuvering  time  and  the  number  of  deployments  per  nominal  coverage.  Assuming 
that  the  vehicle  will  follow  Uuc  paths  of  length 

r,  s  ZJ 

before  turning  onto  adjacent  paths,  there  would  be 

Ml 

r, 

such  turning  maneuvers  during  one  nominal  coverage  of  the  search  area,  and  if  T,  is 
the  time  spent  on  one  turn 


NeJm 
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the  same  scale  factor  multiplying  lengths. 

The  number  of  deployments  is  related  to  the  deployment  time  or  endurance  time  of 
the  vehicle,  denoted  To,  typically  set  by  batteries  or  other  such  factors.  Explicitly 

nNDTD  =  Ts(n)  +  nND(T,  +  T2  +  TM)  +  Tc(n) 

Using  all  of  these  expressions  gives 

„  -  v  [Ts(n)  +  Tc{n)lf^  t  *7*  nML",\ 

T^n)  -  2  "  v 

MEAN  TIME  TO  DETECTION 
Note  first  that 

fxfiwi -*)■-■  -  £  7i^(i ^ 
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These  give 


I' 7b  +  7*  +  71 V Ml  fHLx ^  MLT,  f  Ml  ffHLx y  1  -  <pY 
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PARAMETRIC  STUDIES 


SEARCH  DESIGN  BY  MEAN  TIME  TO  DETECTION 

The  purpose  of  these  parametric  studies  is  twofold,  the  larger  purpose  being  to 
facilitate  the  design  of  a  search  by  quantifying  its  effectiveness  in  terms  of  a  figure  of 
merit,  the  mean  time  to  detection.  Typically  the  target  of  interest  and  the  nature  of  the 
terrain  determine  the  lateral  range  function  and  swath  width,  with  navigation  error 
given  for  a  specific  system  configuration.  Hence,  the  twofold  purpose  of  this  study 
follows  the  nature  of  the  design  of  a  search,  which  consists  of  first  choosing  a  path 
separation  sufficient  to  ensure  an  acceptable  location  probability  for  a  given  lateral 
range  function,  swath  width,  and  navigation  error,  and  then  evaluating  the  effectiveness 
of  this  search  by  the  associated  mean  time  to  detection. 

The  figures  presented  below  are  intended  to  illustrate  the  dependence  of  A ,  A0 , 
location  probability,  and  ultimately  mean  time  to  detection  on  such  system  parameters 
as  swath  width,  navigation  error,  lateral  range  function,  and  path  separation.  The 
validity  of  the  approximation  A0  is  considered  first,  both  in  terms  of  its  closeness  to  A 
and  to  two  associated  location  probabilities,  0  and  .  This  validity  is  important  here 
because  the  approximation  X0  is  used  to  determine  the  desired  path  separation,  as  out¬ 
lined  below.  As  stated  above,  the  figure  of  merit  of  greatest  interest  is  mean  time  to 
detection;  the  relevant  quantity  evaluated  here  is  mean  search  time,  given  as  a  function 
of  swath  width,  navigation  error,  location  probability,  and  lateral  range  function. 


VALIDITY  OF  THE  APPROXIMATION  l0 


Figure  4  shows  the  parameter  A  as  a  function  of  swath  width,  path  separation,  and 
navigation  error  for  a  typical  cookie-cutter  lateral  range  function  3  -  0.9;  the  RMS 
deviation  u  was  taken  to  be  the  RMS  deviation  at  a  point  midway  along  the  nominal 
path,  as  discussed  earlier.  Figure  5  gives  h.  for  the  same  range  of  parameters.  In  both 
figures,  the  quantities  evaluated  are  * 
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with  H,  L  an  order  of  magnitude  greater  than  d,  8,  respectively. 
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5.  X0fX  as  a  function  of  swath  width  for  3  =  0.9. 


Using  Eq.  15  for  X0  and  the  definition  &>  =  l  -  exp 
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for  X0  the  path  separation  d  required  to  ensure  this  location  <j>, 
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gives  an  approximation 
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The  most  important  estimation  of  the  validity  of  the  approximation  X0  is  given  in 
figures  6a~d,  which  show  the  relative  error  between  the  probability  4>  and  a  derived 
location  probability 

*  - 

in  which  the  path  separation  d  is  determined  from  Eq.  16  and  X  is  then  determined 
from  Eq.  14.  These  figures  suggest  using  X0  to  determine  a  path  separation  given  <p , 
(3,  W,  and  ts  is  acceptable  for  a  reasonable  range  of  these  parameters,  where  accept¬ 
ability  is  taken  to  be 

ll0rf-0ll 


s  0.05 


PATH  SEPARATION  AND  MEAN  SEARCH  TIME  PARAMETRICALLY 

In  the  sequence  of  figures  which  follow,  giving  path  separation  (figures  7a-d)  and 
mean  search  time  (figures  8a-d)  as  functions  of  location  probability,  lateral  range 
function,  swath  width,  and  RMS  deviation,  only  those  values  are  shown  which  satisfy 
the  acceptability  criterion 

Mflrf-flll  s  0.05 

0 

The  measure  of  mean  time  to  detection,  which  is  given  in  figures  8a-d,  namely  Vt  is 
derived  from  the  mean  search  time  discussed  above: 

t,  .  -  2*w«*i*(wr,5?r  ■  = 

ov  dv(p 

T  is  a  scale  version  of  Ts- 

T  m 


Note  how  this  scaled  estimate  of  mean  search  time  is  directly  proportional  to  the  path 

ZT  $3 

length  multiplier  a  =—  and  the  value  of  V  actually  plotted  is  More  complete 

studies  of  mean  time  to  detection  for  the  basic  Reber  model  of  the  Advanced 

Unmanned  Search  System  are  included  in  Grace,  1986,  1991. 
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(d)  For  location  probability  ■  0.60. 

Figure  7.  Continued. 
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(d)  For  location  probability  *  0.60. 

Figure  8.  Continued. 
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